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Magnetic helicity in stellar dynamos: new numerical experiments 
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Abstract. The theory of large scale dynamos is reviewed with particular emphasis on the magnetic helicity constraint in 
the presence of closed and open boundaries. In the presence of closed or periodic boundaries, helical dynamos respond to 
the helicity constraint by developing small scale separation in the kinematic regime, and by showing long time scales in 
the nonlinear regime where the scale separation has grown to the maximum possible value. A resistively limited evolution 
towards saturation is also found at intermediate scales before the largest scale of the system is reached. Larger aspect ratios 
can give rise to different structures of the mean field which are obtained at early times, but the final saturation field strength 
is still decreasing with decreasing resistivity. In the presence of shear, cyclic magnetic fields are found whose period is 
increasing with decreasing resistivity, but the saturation energy of the mean field is in strong super-equipartition with the 
turbulent energy. It is shown that artificially induced losses of small scale field of opposite sign of magnetic helicity as the 
large scale field can, at least in principle, accelerate the production of large scale (poloidal) field. Based on mean field models 
with an outer potential field boundary condition in spherical geometry, we verify that the sign of the magnetic helicity flux 
from the large scale field agrees with the sign of a. For solar parameters, typical magnetic helicity fluxes lie around lO*'^ Mx^ 
per cycle. 
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1. Introduction 

The conversion of kinetic into magnetic energy, i.e. the dy- 
namo effect, plays an important role in many astrophysical 
bodies (stars, planets, accretion discs, for example). The gen- 
eration of magnetic fields on scales similar to the scale of 
the turbulence is a rather generic phenomenon that occurs for 
sufficiently large magnetic Reynolds numbers unless certain 
antidynamo theorems apply, which exclude for example two- 
dimensional fields (e.g.. Cowling 1934). 

There are two well-known mechanisms, which allow the 
generation of magnetic fields on scales larger than the eddy 
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scale of the turbulence: the alpha-effect (Steenbeck, Krause 
& Radler 1966) and the inverse cascade of magnetic helic- 
ity in hydromagnetic turbulence (Frisch et al. 1975, Pouquet, 
Frisch, & Leorat 1976). In a way the two may be viewed as 
the same mechanism in that both are driven by helicity. The 
a-effect is however clearly nonlocal in wavenumber space, 
whereas the inverse cascade of magnetic helicity is usually 
understood as a (spectrally) local transport of magnetic en- 
ergy to larger length scales. Both aspects are seen in simula- 
tions of isotropic, non-mirror symmetric turbulence: the non- 
local inverse cascade or a-effect in the early kinematic stage 
and the local inverse cascade at later times, when the field 
has reached saturation at small scales (Brandenburg 2001a, 
hereafter referred to as BOl). 

It is not clear whether either of these two mechanisms is 
actually involved in the generation of large scale magnetic 
fields in astrophysical bodies. Although much of the large 
scale magnetic field in stars and discs is the result of shear 
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and differential rotation, a mechanism to sustain a poloidal 
(cross-stream) magnetic field of sufficiently large scale is 
still needed. The main problem that we shall be concerned 
with in this paper is that of the associated magnetic helic- 
ity production, which may prevent the large scale dynamo 
process from operating on a dynamical time scale. However, 
not many alternative and working dynamo mechanisms have 
been suggested so far In the solar context, the Babcock- 
Leighton mechanism is often discussed as a mechanism that 
operates preferentially in the strong field regime (e.g. Dik- 
pati & Charbonneau 1999). However, Stix (1974) showed 
that the model of Leighton (1969) is formally equivalent 
to the model of Steenbeck & Krause (1969), except that in 
Leighton (1969) the a-effect was nonlinear and mildly sin- 
gular at the poles. Nevertheless, a magnetically driven a ef- 
fect has been suspected to operate in turbulent flows driven by 
the magnetorotational instability (Brandenburg et al. 1995) or 
by a magnetic buoyancy instability (Ferriz-Mas, Schmitt, & 
Schussler 1994, Brandenburg & Schmitt 1998). In the frame- 
work of mean-field dynamo theory, magnetically driven dy- 
namo effects have also been invoked to explain the observed 
increase of stellar cycle frequency with increased stellar ac- 
tivity (Brandenburg, Saar, & Turpin 1998). However, as long 
as these mechanisms lead to an a-effect, they also produce 
magnetic helicity and are hence subject to the same problem 
as before. 

Although kinetic helicity is crucial in the usual expla- 
nation of the a effect, it is not a necessary requirement for 
magnetic field generation, and lack of parity invariance is al- 
ready sufficient (Gilbert, Frisch, & Pouquet 1988). However, 
as in every a effect dynamo, magnetic helicity is produced at 
large scales, which is necessarily a slow process. A mecha- 
nism similar to the ordinary a effect is the incoherent a ef- 
fect (Vishniac & Brandenburg 1997), which works in spite 
of constantly changing sign of kinetic helicity (in space and 
time) provided there is systematic shear and sufficient turbu- 
lent diffusion. This mechanism has been invoked to explain 
the large scale magnetic field found in simulations of accre- 
tion disc turbulence (Brandenburg et al. 1995, Hawley, Gam- 
mie, & Balbus 1996, Stone et al. 1996). In addition, how- 
ever, the large scale field of Brandenburg et al. (1995) shows 
spatio-temporal coherence with field migration away from 
the midplane. This has so far only been possible to explain 
in terms of an afi dynamo with a magnetically driven a ef- 
fect (Brandenburg 1998). In any case, the incoherent a effect, 
which has so far only been verified in one-dimensional mod- 
els, would not lead to the production of net magnetic helicity 
in each hemisphere. 

Yet another mechanism was suggested recently by Vish- 
niac & Cho (2001), which yields a mean electromotive force 
that does not lead to the production of net magnetic helicity, 
but only to a transport of preexisting magnetic helicity. This 
is why this mechanism can, at least in principle, work on a 
fast time scale even when the magnetic Reynolds number is 
large. Finally, we mention a totally different mechanism that 
works on the basis of negative turbulent diffusion (Zheligov- 



sky, Podvigina, & Frisch 2001). This is what produces turbu- 
lence in the Kuramoto-Sivashinsky Equation, which is then 
stabilised by hyperdiffusion (i.e. a fourth derivative term). 
Among all these different dynamo mechanisms, the conven- 
tional a effect and the inverse magnetic cascade are the only 
mechanisms that have been shown numerically to produce 
strong large scale fields under turbulent conditions. These are 
however exactly the mechanisms that suffer from the helicity 
constraint. 

The purpose of the present paper is to assess the signifi- 
cance of the helicity constraint and to present some new nu- 
merical experiments that help understanding how the con- 
straint operates and to discuss possible ways out of the 
dilemma. For orientation and later reference we summarise in 
Fig. |l] various models and the possible involvement of mag- 
netic helicity in them. Dynamos based on the usual a effect 
are listed to the right. They all produce large scale magnetic 
fields that are helical. However, in a periodic or an infinite 
domain, or in the presence of perfectly conducting bound- 
aries, the magnetic helicity is conserved and can only change 
through microscopic resistivity. This would therefore seem to 
be too slow for explaining variations of the mean field on the 
time scale of the 1 1 year solar cycle. Open boundary condi- 
tions may help, although this was not yet possible to demon- 
strate, as will be explained in Sect. Q of this paper 

Given that the issue of magnetic helicity is central to ev- 
erything that follows, we give in Sect. ^ a brief review of 
magnetic helicity conservation, the connection between the 
inverse cascade and the realisability condition, and the sig- 
nificance of gauge-invariant forms of magnetic helicity and 
the surface integrated magnetic helicity flux. Readers famil- 
iar with this may jump directly to Sect. |^, where we discuss 
the issue of magnetic helicity cancellation in kinematic and 
non-kinematic dynamos, or to Sect. ^ where we begin with 
a discussion of the results of BOl and Brandenburg & Dobler 
(2001, hereafter referred to as BDOI), or to Sect. 4.3 where 
new results are presented. 

2. The magnetic helicity constraint 

Magnetic helicity evolution does not give any additional in- 
formation beyond that akeady contained in the induction 
equation governing the evolution of the magnetic field it- 
self. Nevertheless, the concept of magnetic helicity proves 
extremely useful because magnetic helicity is a conserved 
quantity in ideal MHD, and almost conserved even in slightly 
non-ideal conditions. The evolution equation for the mag- 
netic helicity can then be used to extract information that 
can be easily understood. At a first glance, however, mag- 
netic helicity appears counterintuitive, because it involves the 
magnetic vector potential. A, which is not itself a physical 
quantity, as it is not invariant under the gauge transformation 
A ^ A — VV'. Only the magnetic field, B = V x A, which 
lacks the irrotational 'information' from A, is independent 
of t/) and hence physically meaningful. Therefore A • B is 
not a physically meaningful quantity either, because a gauge 



2 



Astron. Nachr./AN XXX (200X) X 



Large scale dynamos 

(e.g. in the sun) 



Is significant net magnetic helicity 
produced in eacli hemispliere? 



Using Eqs (^) and one can derive an equation for the 
gauge-dependent magnetic helicity density A • B, 

d 

— (A- B) = -2E- B- V • (E X A + ^-B). (5) 

Using Ohm's law, rjfioJ = E + u x B, we have 2E • B = 
2ry/io J • B, where rj is the microscopic magnetic diffusivity, 
fiQ the magnetic permeability in vacuum, and J = V x B//^o 
the electric current density. 
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Fig. 1. Sketch illustrating the different ways astrophysical dy- 
namos may be able to circumvent the magnetic helicity prob- 
lem. 

transformation, A — > A — 'Vtp changes A • B. For periodic 
or perfectly conducting boundary conditions, or in infinite 
domains, the integral H ~ J A • B dV is however gauge- 
invariant, because in 



IV JdV Jv 

the surface integral on the right hand side vanishes, and the 
last term vanishes as well, because V • B = 0. 

The evolution of the magnetic field is governed by the 
induction equation 
dB „ 

^ = X E, (2) 
ot 

where E is the electric field. This equation can be integrated 

to yield 

(9A 

— = -E - V$, (3) 
ot 

where $ is the scalar ('electrostatic') potential, which is also 
refeiTed to as the gauge potential, because it can be chosen 
arbitrarily without affecting B. Common choices are V^"I> = 

— V • E, which preserves the Coulomb gauge V • A = 0, and 
$ = 0, which is convenient in MHD for numerical purposes. 
If 77 = const, another convenient gauge is $ = —77V • A, 
which results in 

dA 

— — = u X B + nV^A (77 = const). (4) 

ot 



2.1. Magnetic helicity conservation with periodic 
boundaries 

Consider first the case of a periodic domain, so the divergence 
term in Eq. (||) vanishes after integration over the full volume, 
and thus 

A(A.B) = -27;/io(J-B), (6) 



where angular brackets denote volume averages, so J A • 
B dV = ( A • B) V, where V is the (constant) volume of the 
domain under consideration. Note that the terms in Eq. (^ 
are gauge-independent [see Eq. (^]. More important, how- 
ever, is the fact that the rate of change of magnetic helicity is 
proportional to the microscopic magnetic diffusivity rj. This 
alone is not sufficient to conclude that the magnetic helic- 
ity will not change in the limit 7; 0, because the current 
helicity, (J • B), may still become large. A similar effect is 
encountered in the case of ohmic dissipation of magnetic en- 
ergy which proceeds at the rate (Jjouic = '7/^0 (J^)- In the 
statistically steady state (or on time averaging), the rate of 
ohmic dissipation must be balanced by the work done against 
the Lorentz force. Assuming that this work term is indepen- 
dent of the value of 7; (e.g. Galsgaard & Nordlund 1996) we 
conclude that the root-mean-square current density increases 
with decreasing rj like 

-^rms OC 77^1/2 as 7/ 0, (7) 

whilst the rms magnetic field strength, Brms, is essentially in- 
dependent of 77. This, however, implies that the rate of mag- 
netic helicity dissipation decreases with rj like 

Qff = 277/io(J ■ B)I/ cx 77+1/2 ^ as 77 ^ 0. (8) 

Thus, under many astrophysical conditions where the mag- 
netic Reynolds number is large (77 small), the magnetic helic- 
ity H, as governed by Eq. (||), is almost independent of time. 
This motivates the search for dynamo mechanisms indepen- 
dent of H. On the other hand, it is conceivable that magnetic 
helicity of one sign can change on a dynamical time scale, 
i.e. faster than in Eq. if magnetic helicity of the other sign 
is removed locally by advection, for example. This is essen- 
tially what the mechanism of Vishniac & Cho (2001) is based 
upon. So far, however, numerical attempts to demonstrate the 
operation of this mechanism have failed (Arit & Branden- 
burg 2001). Yet another possibility is that magnetic helicity 
of one sign is generated at large scales, and is compensated 
by magnetic helicity of the other sign at small scales, as in 
the kinematic a effect. However, as we shall see in Sect. ^ 
this does not happen efficiently when nonlinear effects of the 
Lorentz force come into play. 
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2.2. Realisability condition and connection with inverse 
cascade 

Magnetic helicity is important for large scale field genera- 
tion because it is a conserved quantity which cannot easily 
cascade forward to smaller scale. We shall demonstrate this 
here for the case where the magnetic field is fully helical. 
The existence of an upper bound for the magnetic helicity 
is easily seen by decomposing the Fourier transformed mag- 
netic vector potential, Aj^^, into a longitudinal component, 
ft.ll, and eigenfunctions of the curl operator (also called 
Chandrasekhar-Kendall functions). 



with 



fly IXi IX, fly IXi IXi 



Wxhl^±khl, k^\k\, 
and 

{h+* ■h+) = {h7; ■hj^) = {hi-hi) = i, 



(9) 



(10) 



(11) 



where asterisks denote the complex conjugate. The longitu- 
dinal part is parallel to k and vanishes after taking the 
curl to calculate the magnetic field. In the Coulomb gauge, 
V • A = 0, the longitudinal component vanishes altogether 
(apart from an uninteresting constant vector Aq). Since 
and ftll are all orthogonal to each other, it is clear that, for a 
given magnetic field, the Coulomb gauge minimises (A^) — 
or, if Ao = (A) ^ 0, the variance {{A-{A)f). 

The (complex) coefficients a^(i) depend on k and t, 
while the eigenfunctions h^, which form an orthonormal set, 
depend only on k and are given by[| 

^± ^ 1 fc X (fc X e) Tifc(fc X e) 
fe"V2 k^y^T~{k~ejyk^ ' 

where e is an arbitrary unit vector not parallel to fc. With 
this definition we can write the magnetic helicity and energy 
spectra in the form 



k{\a 



+ |2 



)V, 



)V, 



(13) 
(14) 



where V is the volume of integration. (Here and elsewhere 
the factor /Iq ^ is ignored in the definition of the magnetic 
energy.) These spectra are normalised such that 



Hkdk^ (A-B)y = iJ, 



Mi.dfc 



'1t32 



B^)V=M, 



(15) 



(16) 



where H and M are magnetic helicity and magnetic energy, 
respectively. From Eqs ( |l3| ) and ( [Til ) one sees immediately 
that 



^k\Hk\ < Mfc, 



(17) 



The forcing functions used in BOl were proportional to hj^. We 
note that, by mistake, the square root was missing in the normalisa- 
tion factor 1 / \/2 in that paper 



which is also known as the realisability condition. A fully 
helical field has therefore Alk — ±\kHk. 

For the following it is convenient to write the magnetic 
helicity spectrum in the form 



Hk = {Ak-Bk)V/8k, 



(18) 



where the subscript k (which is here a scalar!) indicates 
Fourier filtering to only those wave vectors fc that lie in the 
shell 

k~ 5k/2>\k\<k + 5k/2 (fc-shell). (19) 
Note that A^ and are in real (configuration) space, so 



V X Afc. (By contrast, Aj^ and Bj^ = ifc x A 



fe 



are quantities in Fourier space, as indicated by the bold face 
subscript fc, which is here a vector) In Eq. ( [l8| ) the average 
is taken over all points in x space. Note that in practice the 
helicity is more easily calculated directly as 



(20) 



''^^^Ishel/^^-'''^ + ^^-^^^^'''''^' 
where AVlj^ is the solid angle element in Fourier space. The 
result is however the same. Likewise, the magnetic energy 
spectrum can be written as 



Mk = h{'Ql)Vl5k, 



or as 



1 



B 



fc 



B^ fc2 dD.k 



(21) 



(22) 



2 Jfc-shell 

We recall that for a periodic domain H is gauge invariant. 
Since its spectrum can be written as an integral over all space, 
see Eq. ([T^, Hk is - like H - also gauge invariant. 

The occurrence of an inverse cascade can be understood 
as the result of two waves (wavenumbers p and q) interact- 
ing with each other to produce a wave of wavenumber fc. The 
following argument is due to Frisch et al. (1975). Assum- 
ing that during this process magnetic energy is conserved to- 
gether with magnetic helicity, we have 

(23) 



Mp + Mg = Mk 



\H,\ - \Hk\ 



(24) 



where we are assuming that only helicity of one sign is in- 
volved. Suppose the initial field is fully helical and has the 
same sign of magnetic helicity at all scales, then we have 

2Mp=p\Hp\ and 2M, = q\H,\, (25) 



and so Eq. (23 ) yields 

p\Hp\ + q\H,\ = 2Mk > k\Hkl (26) 

where the last inequality is just the realisability condition ( |T7[ ) 
applied to the target wavenumber fc after the interaction. Us- 
ing Eq. ( |24| ) we have 

p\Hp\+q\Hg\>k{\Hp\ + \Hg\). (27) 

In other words, the target wave vector fc after the interaction 
of wavenumbers p and q satisfies 

p\Hp\+q\H,\ 



k < 



(28) 
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The expression on the right hand side of Eq. ( p8| ) is a 
weighted mean of p and q and thus satisfies 

min(p, q) < ^l^^'+flf^' < max(p, q), (29) 

and therefore 

k < max(p, g). (30) 

In the special case where p = q,we have k < p = q, so the 
target wavenumber after interaction is always less or equal to 
the initial wavenumbers. In other words, wave interactions 
tend to transfer magnetic energy to smaller wavenumbers, 
i.e. to larger scale. This corresponds to an inverse cascade. 
The realisability condition, ^k\Hk\ < Affc, was perhaps the 
most important assumption in this argument. Another impor- 
tant assumption that we made in the beg inning was that the 
initial field be fully helical; see Sect. 4.5 for the case of frac- 
tional helicity. 

We note that in hydrodynamics, without magnetic fields, 
the kinetic helicity spectrum, Fk ~ {u>k ■ Uk)V/6k, where 
u) = V X u is the vorticity, is also conserved by the nonUnear 
terms, but now the realisability condition is 



(31) 



so it is not fc, but the inverse of k, that enters the inequality. 
Here, Ek ^ \ {^)^ related to the kinetic energy spectrum 
(i.e. without density), and hence one has 



< 



if; 



Fa 



or 



k ^ < max(p ^ ,q ^ 



(32) 



(33) 



Thus k > mm{p, q), and in the special case p — q we have 
k > p — q, so kinetic energy is transferred preferentially 
to larger wavenumbers, so there is no inverse cascade in this 
case. 



2.3. Spectra of right and left handed components 

For further reference we now define power spectra of those 
components of the field that are either right or left handed, 
i.e. 



= ±k\a±\^V, Af± = ^k^\a±\'^V. 



(34) 



Thus, we have Hk = + and Mk = + . Note 
that and can be calculated without explicit decom- 
position into right and left handed field components using 

= \{Hk ± 2k-^Mk), M± = ± \kHk). (35) 

This method is significantly simpler than invoking explicitly 
the decomposition in terms of a^^^- 



2.4. Open boundaries and the importance of a 
gauge-independent magnetic heUcity 

Let us now discuss the case where magnetic helicity is al- 
lowed to flow through boundaries. The first problem we en- 
counter in connection with open boundaries is that the mag- 
netic helicity is no longer gauge invariant. There are how- 
ever known procedures that allow H to be defined in a way 
that is independent of gauge. Before we discuss this in more 
detail we first want to demonstrate how the use of a gauge- 
dependent magnetic helicity can lead to undesired artifacts. 

For illustrative purposes we consider the case of a one- 
dimensional mean-field dynamo governed by the equation 



(36) 



where ?7x = r^t + ?y is the sum of turbulent and microscopic 
magnetic diffusivity, and B = {B^^By^O) — B(z) is the 
magnetic field which depends only on z, and B = V x A 
as well as J = V x B//io- Here we have used the gauge 
$ = 0. The overbars denote horizontal (x, y) averages and 
the corresponding mean quantities like B are therefore inde- 
pendent of X and y. As boundary conditions we adopt either 
dAx/dz = dAy/dz = (vacuum boundary condition) or, 
for comparison, = Ay = (perfect conductor condition) 
on z = 0, L^. We adopt a and rjt quenching, i.e. 

"0 „ ?/to 



1 



aB^'/Bl 



Vt 



cq 



1 



(37) 



cq 

(In the case considered here we use a = 1.) 

The result of a numerical integration of Eq. ( ^6| ) is shown 
in Fig. ^ Evidently, the use of a gauge-dependent magnetic 
helicity implies spurious contributions that prevent (A • B) 
from becoming time-independent even though the magnetic 
energy has long reached a steady state. The reason for this 
is that for the steady state solution the right hand side of 
Eq. (^6|) is a constant which is different from zero. As an 
example, consider the eigenfunction of the linearised equa- 
tions: B = (1 + coskiz,smkiz,0)BQ, where ki = 1 for 
Lz = 2tt. In the marginally excited state, we have a ~ r^xfci, 
sodAx/dt = aBo{t) ^ 0. Because of this complication, it is 
important to define magnetic helicity in a gauge-independent 
way. BDOl defined such a gauge-independent magnetic he- 
licity by linearly extending the horizontally averaged mag- 
netic vector potential to a periodic field which then results 
in a gauge-invariant magnetic helicity as before; see Eq. ([l|). 
This can be written in compact form as 



BD 



(A-B>y + i-(Ai xA2)V, 



(38) 



where Ai and A2 are the values of A on the lower and up- 
per boundaries, and angular brackets denote full volume av- 
erages, as opposed to the overbars, which denote only hori- 
zontal averages.^ 

^ We use this opportunity to point out a sign error in Eq. (9) of 
BDOl, where it should read Ao — Ao — Vx x (i^z). The contri- 
butions from the mean field, that are central both here and in BDOl, 
remain however unaffected. 
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Fig. 2. First panel: gauge-dependent magnetic helicity for an 
dynamo with vacuum boundary conditions. Second panel: 
gauge-independent magnetic helicities H^d (solid line) and 
^^Cou (dashed line), compared with twice the magnetic en- 
ergy, 2AI, for the same model (dotted line). Third panel: same 
as second panel, but with perfect conductor boundary condi- 
tions. Here gauge-independent magnetic helicities agree with 
the ordinary expression, J A • B dy. Again, the dotted line 
shows 2M. 

We note that Hbu can be written in the form proposed by 
Finn &Antonsen (1985), 

irBD(B) = ((A + Ap).(B-Bp))F (39) 
where Bp is the non-helical reference field 
Bp = (B), _ (40) 

which is spatially constant. A corresponding Ap that satisfies 
Bp = V X Ap is 

Ap = — z X Bp + const, (41) 



which, apart from some constant, corresponds to a linear in- 
terpolation of A between zi and Z2- However, the constant 
in Eq. ( |4l| ) drops out in the term Ai x A2 of Eq. (Esh. It is 
also easy to see that the constant drops out in Eq. (39), be- 
cause (B — Bp) = 0. It is precisely this fact that makes ( ^9| ) 
gauge-invariant. 

In the following we compare with the magnetic helicity 
A), which minimises (A ), 



in the gauge A i-^ A 
i/cou-(A.B)l/-(A: 



(B)V. 



(42) 



This magnetic helicity is actually also gauge-invariant, but it 
cannot be written in the form (P9|), because it is not based on 
a reference field. Indeed, the only value of Bp that makes the 
expression ( ^9| ) gauge-invariant is Bp = (B), but this does 
already corresponds to Hbu- 

We associate -ffcou with the Coulomb gauge because it 
is based on a gauge which minimises (A^), which is char- 



acteristic of the Coulomb gauge; see Sect. |2.2[ However, in 
the present geometry, for horizontally averaged fields, and 
zero net magnetic flux, a canonical gauge transform A 1— s- 
A — zdtp I dz does not change the large scale magnetic he- 
licity (A • B), since Bz ~ 0. In this sense, (A • B)t/ in any 
arbitrary gauge represents the magnetic helicity in Coulomb 
gauge. There is, however, the extra freedom to add a constant 
vector — Aq to the magnetic vector potential, A A — Aq. 
It is easy to see that (A ) is minimised if Aq = (A), and that 
the definition ( p^ always yields the magnetic helicity which 
corresponds to this gauge. 

In the second panel of Fig. || we have plotted both iJeo as 
well as Hcon- Note that |i?cou| < \Hbi} \ and also |iJcou| < 
2Af/fcniin. If the domain was periodic, and since the compu- 
tational domain is — tt < a; < tt, we would have fcmin — 1. 
However, for vacuum boundary conditions one could still ac- 
commodate a wave with k^in = 1/2. Then also |i?BD| is 
always less than 2Af/fcmin. 

In the more general case of three-dimensional fields we 
adopt the gauge-independent magnetic helicity of Berger & 
Field (1984). Their definition, however, leaves the contribu- 
tion from the average of A undetermined and can therefore 
only be applied to the deviations from the mean field. 

2.5. Resistively limited growth and open boundaries 

A plausible interpretation of the resistively limited growth is 
the following (Blackman & Field 2000, Ji & Prager 2001). 
To satisfy magnetic helicity conservation, the growth of the 
large scale field has to proceed with very little (or very slow) 
changes of net magnetic helicity. This is essentially the result 
of BOl (see also Brandenburg & Subramanian 2000, hereafter 
BSOO). The underlying assumption was that magnetic helic- 
ity of the sign opposite to that of the mean field is lost resis- 
tively. As a working hypothesis we assume that such losses 
occur preferentially at small scales where also Ohmic losses 
are largest. 

In order to speed up this process it might help to get rid of 
magnetic helicity with sign opposite to that of the large scale 
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field. One possibility would be via explicit losses of field with 
opposite magnetic helicity through boundaries (i.e. the outer 
surface or the equatorial plane). This possibility was first dis- 
cussed by Blackman & Field (2000) and more recently by 
Kleeorin et al. (2000) and Ji & Prager (2001). 

The results obtained so far seem to indicate that this does 
not happen easily (BDOl). Instead, since most of the mag- 
netic helicity is already in large scales, most of the losses are 
in large scales, too. Thus, the idea of losses at small scales, 
where the sign of magnetic helicity is opposite to that at large 
scale, has obvious difficulties. In Sect. 4.1 we address the 
question of observational evidence for this. In Sect. ^ we dis- 
cuss quantitatively the effect of losses on large scales. 

3. Magnetic helicity cancellation 

3.1. Magnetic helicity in ABC-flow dynamos 

Before we start discussing the evolution of magnetic helic- 
ity in turbulent dynamos we consider first the case of steady 
kinematic ABC flows, where the velocity field is given by 



u 



As'mkiz + C coskiy 
B sin kix + A cos kiz 
. C sin kiu + B cos kix ^ 



(43) 



with A = B = C = 1 and ki = 1 in the domain 
—77 < x,y,z < IT. This class of dynamo models is of in- 
terest because here the net magnetic helicity can be small, 
even though the flow itself is fully helical (Gilbert 2002, and 
references therein). The field generated by such an ABC flow 
dynamo does have magnetic helicity, but it is comparatively 
weak in the sense that the realisability condition ( |l7[ ) is far 
from being saturated. This is seen in Fig. |], where we com- 
pare Mfc with i k I Hk I . In addition, the spectral helicity is al- 
ternating in sign, often from one wavenumber to the next; 
see the second panel of Fig. |] The magnetic helicity den- 
sity is also alternating in space, with a smoothly varying, 
mostly positive component, interspersed by negative, highly 
localised peaks as shown in Fig. ^ Moreover, as oo, 
the net magnetic helicity goes asymptotically to zero; see 
Fig. H. This result is typical of fast dynamos (Hughes, Cat- 
taneo, & Kim 1996, Gilbert 2002). 

We note in passing that the magnetic power spectrum in- 
creases approximately like k^^^ (see the left panel of Fig. 
which is similar to the results from a number of other dy- 
namos as long as they are in the kinematic regime; see, e.g.. 
Fig. 22 of BOl for a large magnetic Prandtl number calcu- 
lation and Fig. 17 of Brandenburg et al. (1996) for convec- 
tive dynamos. Thus, most of the spectral magnetic energy is 
in small scales, which is quite different from the helical tur- 
bulence simulations of BOl. In turbulent dynamos this spec- 
trum could be explained by the analogy between magnetic 
field and vorticity (Batchelor 1950). For a Kolmogorov k^^/^ 
spectrum, the vorticity spectrum has a k^/'^ inertial range. In 
the present case, however, this argument does not apply be- 
cause the flow is stationary, so there is no cascade in velocity. 




5 10 15 20 
k 

Fig. 3. Power spectrum, Af^, of the magnetic field for an 
ABC flow dynamo with A = B = C=Uo^l and 
fc = 1. Note that most of the power is in the small scales with 
Mk ^ k^/^. The normalised helicity spectrum is also plotted. 
Note that for medium and large values of k the spectral mag- 
netic helicity is very small; \k\Hk\ ^ Mk- The simulation 
was done with 240'^ mesh points and a magnetic diffusivity 
?7 = 2 X 10-4. 




Fig. 4. Cross-section of the magnetic helicity density for the 
ABC flow dynamo from Fig. |[ Note the intermittent nature 
of AB. 



Thus, the origin of the k^^^ spectrum for these types of kine- 
matic dynamos remains unclear 

One might wonder whether nearly perfect magnetic helic- 
ity cancellation can also occur in the sun. Before we can ad- 
dress this question we need to understand the circumstances 
under which this type of magnetic helicity cancellation can 
occur It appears that nearly perfect magnetic helicity cancel- 
lation is an artefact of the assumption of linearity and of a 
small computational domain. Increasing its size allows larger 
scale fields to grow even though the field has already satu- 
rated on the scale of the ABC flow. This is shown in Fig. ^ 
where we plot the magnetic energy of the mean magnetic 
field, B, for a 'nonlinearised' ABC flow with A — and 
B ~ C — 1/(1 + Af /Mq), where M is the magnetic energy 
and Afo = 1- We chose a domain of size 2tt x 2tt x Stt, 
which is increased in the z direction by a factor 4. In the 
kinematic phase the high wavenumber modes grow fastest. 
Once the magnetic energy reaches saturation, the velocity is 
decreased, so the high wavenumber modes stop growing, but 
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Fig. 5. Magnetic helicity, (A • B), nondimensionalised by 
(B'^)/A;i, as a function of i?in = [/q/ (f?^! ) for a model with 
A — B — C=Uq~1. Here, fci = 1 is the smallest 
wavenumber in the box, which has the size {27t)^. Note that 
A:i(A • B)/(B^) decreases with increasing like i?n/^^- 
Thus, in the high magnetic Reynolds number limit 1/?/ ^ oo, 
the magnetic helicity of the ABC flow dynamo vanishes. 

the remaining lower wavenumber modes are still excited. Al- 
though they grew slower during the earlier kinematic stage, 
they are now still excited and saturate at a much later time. 
This process continues until the mode with the geometrically 
largest possible scale has reached saturation, provided that 
mode remains still excited. Once a large scale field of this 
form has developed, cancellation of any form has become 
small. This is indicated by the fact, that the magnetic helic- 
ity is nearly maximum, i.e. H/ (2/ioA/) w k^^^ with km be- 
ing the momentary wavenumber of the large scale magnetic 
field, which is also indicated in Fig. ^. Thus, for this type of 
dynamo, small scale magnetic fields are highly intermittent 
and can grow fast, while large-scale fields grow once again 
on slow time scales only. It is this type of large scale field 
generation that is often thought to be responsible for generat- 
ing the large scale field of the sun (see Parker 1979, Krause & 
Radler 1980, or Zeldovich et al. 1983, for conventional mod- 
els of the solar dynamo), although the saturation process is 
still resistively slow. 

3.2. Magnetic helicity cancellation in turbulent dynamos 

Magnetic helicity cancellation seems to be a generic phe- 
nomenon of all kinematic dynamos, including the turbulent 
dynamos in cartesian domains that have been studied recently 
in BOl. Even though the initial field is fully helical, the signs 
of magnetic helicity are different at large and small scales 
and can therefore cancel, at least in principle. In practice 
there remains always some finite net magnetic helicity, but 
as the magnetic Reynolds number increases, the total mag- 
netic helicity decreases, exa ctly like for the ABC-flow dy- 
namos studied in Sect. 



3.1 



Examples of magnetic power 
spectra calculated separately for the right and left handed 
parts of the field are shown in Fig. |7] for Runs A, B, and C of 



20 



15 - 



fc = 0.5 



10 



k=1.0 
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t 

Fig. 6. Evolution of the magnetic energy M — (B^/2/io), 
for a 'nonlinearised' ABC flow with A = Q, and B = C = 
1/(1 + M/Mq) scales with the magnetic energy, M. The 
box has a side length of Stt in the z direction. In the kine- 
matic stage (t < 500) modes with fc > 1 grow fastest. Af- 
ter the small scale modes saturate, the energy of the larger 
scales continues to grow until the large scale reaches satura- 
tion (e.g., k = 0.5 saturates after t = 2000). The field contin- 
ues to grow even further, but very slowly. Since the box has a 
side length of Stt, the minimum wavenumber is fc = 0.25. 



Brandenburg & Sarson (2002, hereafter referred to as BS02). 
These are models with helical forcing around a wavenumber 
fcf, like in BOl. The basic parameters of the runs presented 
in BS02 are summarised in Table |l|. Run A has ordinary 
magnetic diffusivity with rj — 10^*, whilst Runs B and C 
have second order hyperdiffusivity, so rjW^A is replaced by 
— 772V^A in the induction equation (Q), and 772 = 3 x 10^® 
and 10~^ for Runs B and C, respectively. [Generally, for 
models with nth order hyperdiffusivity, one would have a dif- 
fusion term (— 1)"+^77„V^"A.] In all three cases the forcing 
of the flow is at wavenumber fcf = 27 with the same am- 
plitude (/o = 0.1), as in BOl. At the forcing scale. Run B 
is by a factor r/ifcj /(7/2^f ) — 0.22 times less resistive than 
Run A. Run C is another three times less resistive. During 
the kinematic phase shown, the spectral power increases at 
all wavenumbers at a rate proportional to e'''*, where j — 2X 
is the kinematic growth rate of the magnetic energy and A the 
growth rate of the rms magnetic field. 

As the magnetic Reynolds number increases, the mag- 
netic helicity normalised by the magnetic energy, —H/M, 
decreases, as shown in Fig. ^ for Runs A, B, and C of BS02. 
The helicity of the forcing and the kinetic helicity of the flow 
are positive, so the magnetic helicity at small scales is also 
positive. This results in large scale magnetic helicity that is 
negative. Therefore, also the total magnetic helicity is nega- 
tive. 

The asymptotic decrease of magnetic helicity for progres- 
sively smaller magnetic diffusivity (Fig. ^) can be quantita- 
tively described by Berger's (1984) magnetic helicity con- 
straint, which we modify here for the initial kinematic evolu- 
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Fig. 7. Power spectra separately for the right and left handed 
parts of the field, for Runs A, B, and C of BS02. Different 
times have been collapsed onto the same graph by scaling the 
spectra with e^''*, where 7 is the corresponding growth rate 
of the magnetic energy for each run. Note that as the mag- 
netic diffusivity is reduced, the relative excess of over 
decreases, as expected for asymptotic magnetic helicity 
conservation. 

tion when the magnetic energy increases exponentially at the 
rate 7 — 2A. We consider first the case of ordinary magnetic 
diffusivity and turn then to the case of hyperdiffusivity. For a 
closed or periodic domain, the modulus of the time derivative 
of the magnetic hehcity can be bounded from above by 

1 AH 

Mo 



(44) 



= 27y|(J-B)|y < 2r?(j2)i/2(B2)i/2^ 
- 2^^/2Qj/^j2M)V2, 
— (77/ioJ^) is the rate of Joule dissipation 



where Qjouio 

In the kinematic regime, \H\ = j\H\ and M ~ 7^1/ = 
W — Qjouic = eQjouic, where W is the work done on the 
magnetic field by the velocity field. In the dynamo case, i.e. 
if 7 > 0, we have e > 0. Thus 

— \H\ < 277^/2(^^^/^)1/2 (2^/) 1/2 
Mo 
or 



(27y7/e)i/2(2Af)(45) 



\H\/{2^loM) < e-i/2(2r;/7)i/2 = a4 



(46) 



where ^skin = (277/7) is the skin depth associated with 



the growth rate 7 and a 



-1/2 



is a parameter describing the 



fractional excess of magnetic energy gain over dissipation. 

In the case of hyperdiffusivity, which is often used for 
numerical reasons, we have 

'B2^ 



; Joule 



Vn 



|fc|2ni^rf3^ 



2n 



V, (47) 



Mo \ Mo 

where fcf cff ~ fcf is an effective wavenumber defined such 
that the second equality in Eq. (^8|) below is satisfied exactly 
and in Eq. only approximately. The rate of change of 
magnetic helicity is then bounded by 

^\H\ < 277„/|fe|2»-i|i3^|2/;iod'fc 

Mo , (48) 



or 

-1^1 
Mo 



= (277„A:f2;-2)i/2(^/,)i/2(2M). 



(49) 
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Fig. 8. Magnetic helicity, normalised by the magnetic energy, 
for Runs A, B, and C of BS02. In Run A, 7; = 10"'', whereas 
in Runs B and C second order hyperresistivity is used, with 
772 = 3 X 10^^ and 10"^, respectively. Note that the smaller 
the magnetic diffusivity (at large and intermediate scales), the 
slower is the saturation of magnetic helicity. 



So, again, 

\H\/i2M) < a4kin, 



(50) 



but now £skin ~ (2'7nfct off ll) involves the value of 
fcf^cff- Like in the case with ordinary magnetic diffusivity we 
have a — e"^/^. From the values of |i7|/(2//o^^) and ^skin 
given m Table [l| we see that this constraint is indeed well sat- 
isfied during the kinematic growth phase. 

3.3. Magnetic helicity conversion to large scales 

At least initially, i.e. when the field is still growing exponen- 
tially, there is a strong conversion of positive magnetic he- 
licity at the driving scale to negative magnetic helicity at a 
certain larger scale. This can be described by the magnetic 
helicity equation, written separately for positive and negative 
parts, 

H±=±2S -2f]k'iH±, (51) 

where H± ^ J H^dk, with k^ = J /fc"ff^dfc/ / H^dk, 
and S describes the transfer of magnetic helicity from small 
to large scales. In runs with second order hyperresistivity, the 
term rjkj. has to be replaced by 772 fc± . The evolution of the net 
magnetic helicity, H^-\-H~ , is of course unaffected by S. We 
note, however, that the spectrum of the magnetic helicity does 
have non-diffusive source and sink terms such that its integral 
over all wavenumbers vanishes (see BOl for detailed plots). 
The term ±5 in Eq. ( |5l| ) is related to the a effect (Blackman 
& Field 2000) or, more precisely, to the turbulent electromo- 
tive force. In the present case, where turbulent diffusion also 
enters, we have therefore (cf. BOl) 



S^{aB -77tMoJ-B)y. 



(52) 



Since the other two terms in Eq. ( pl| ) are known, we can di- 
rectly work out S from this equation. The results for Runs A, 
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Table 1. Summary of runs, some of which were already presented in BS02. The runs with n = 2 have hyperdiffusivity, so 
the effective diffusivity at the forcing wavenumber, fcf, and at the Nyquist wavenumber, fcNy = tt/Sx — N/2, are different, 
is the total number of mesh points, and Sx is the mesh spacing. The parameter £skin gives an approximate upper bound 

for \H\/{2^ioM). 



Run 


iV 


n 






2 


'/II i^Ny 


2 


7 


ki 




\H\/{2tioM) 


^skin 


S/{2M)\t= 


A 


120^ 


1 




10"* 




10"" 




0.047 


27 


1.1 X 10"^ 


0.035 


0.065 


0.00075 


B 


120^ 


2 


3 X 10"* 


2 X IQ- 


-5 


10"* 




0.070 


27 


7.3 X 10"* 


0.018 


0.025 


0.00035 


C 


120^ 


2 


10"* 


7 X IQ- 


-6 


4 X IQ- 


-5 


0.082 


27 


3.6 X 10"* 


0.005 


0.013 


0.00020 


D 


30^ 


2 


10"^ 


9 X IQ- 


-4 


2 X 10" 


-2 


0.078 


3 




0.08 


0.15 




E 


60^ 


2 


lO"'^ 


2 X 10~ 


-4 


9 X IQ- 


-3 


0.12 


5 




0.040 


0.065 




F 


120^ 


2 


3 X 10"* 
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-7 


10"* 
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Fig. 9. Relative magnetic helicity transfer from small to large 
scales for runs A (most resistive), B and C (least resistive) of 
BS02, as a function of time. Note that for smaller magnetic 
diffusivity, the ratio S/ {2M) decreases and hence the conver- 
sion to large scales becomes less efficient. 

B, and C are plotted in Fig. I The quantity S'/(2M) IS a mea- 
sure of the rate of helicity transfer, i.e. the difference between 
a-effect and turbulent diffusion. Although the ratio S/ (2Af ) 
has not yet settled to a stationary value, we compare in Ta- 
ble |l] the values for three different runs (at t = 1000). 

Figure ^ shows that S/ {2AI) approaches a finite value in 
the nonlinear regime, but its value becomes smaller as the 
magnetic diffusivity is decreased. Thus, even though strong 
final field strengths are possible, this can still be achieved 
with a residual a effect that decreases with decreasing mag- 
netic diffusivity. This is consistent with earlier results of Cat- 
taneo & Vainshtein (1991), Vainshtein & Cattaneo (1992), 
and Cattaneo & Hughes (1996) that a and i-jt are quenched 
in a i?in -dependent fashion. Field & Blackman (2002) have 
pointed out, however, that this result can also be explained by 
a dynamically quenched a where a and rjt are not necessarily 
quenched in a i?ni -dependent fashion. 

One may wonder whether the i?m-dependent slow-down 
during the saturation phase of the large scale field (as seen in 
Fig. applies also to intermediate scales. In Fig. |l^ we show 
how the magnetic power spectrum shows a bump at interme- 
diate wavenumbers (i.e. the wide maximum at wavenumbers 




1 10 

Fig. 10. Evolution of magnetic energy spectra in the hyper- 
diffusive Run B in equidistant time intervals from t — (low- 
est curve), t — 4:0 (next one higher up), until t — 2080 (peak- 
ing at the very top left). The dotted lines give the kinetic en- 
ergy. 



k < 20fci) travelling to the left toward smaller wavenum- 
bers. In order to check whether or not the speed of the spec- 
tral bumps in the different runs shown in Figs ^ and |l^ 
also decreases with decreasing magnetic diffusivity we plot 
in Fig. |ll] the position of the secondary bump in the power 
spectrum as a function of time. It is evident that the slope, 
which signifies the speed of the bump, decreases with de- 
creasing magnetic diffusivity. A reasonable fit to this migra- 
tion is given by 



(53) 



where the parameter atrav characterises the speed at which 
this secondary bump travels. The values obtained from fits to 
various runs are also listed in Table |l} Following arguments 
given by Pouquet et al. (1976), atrav obtained in this way is 
actually a measure of the a effect. The decrease of atrav with 
decreasing magnetic diffusivity adds to the evidence that the 
a effect is quenched in a i?m dependent fashion. 

The dependence of atrav ™d S/{2M), both indica- 
tive of inverse energy transfer, seems to be in conflict with 
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Fig. 11. Evolution of the position of the secondary bump in 
the power spectrum. 
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Fig. 12. Evolution of magnetic energy in runs with mag- 
netic driving. Forcing occurs in the induction equation around 
wavenumber fcf = 5. Runs with hyperdiffusivity are denoted 
by r\2 and those with ordinary diffusivity by 771 . 



that the magnetic energy evolves to larger values as the mag- 
netic diffusivity is decreased. This is because the injection 
of magnetic energy is primarily balanced by magnetic dif- 
fusion, which becomes weaker as the magnetic diffusivity is 
decreased. In the runs with second order hyperdiffusivity (de- 
noted by r\2) the magnetic energy evolves to even larger val- 
ues. What is more important, however, is the time scale on 
which saturation is achieved. This time scale is the resistive 
one, suggesting again that the transport of magnetic energy 
across the spectrum from small to large scales is a resistively 
limited process. 



4. Evolution of the large scale magnetic helicity 



the usual idea of an inverse cascade where such transfers hap- 
pen on a dynamical time scale (e.g., Verma 2001). The only 
reasonable conclusion from this seems to be that the (nonlo- 
cal!) inverse transfer discussed here is in fact distinct from the 
usual inverse cascade. In other words, the present simulations 
do not support the concept of an inverse cascade in the usual 
(local) sense. 

3.4. Magnetically driven turbulence 

In the cases considered so far, the turbulence is driven by a 
forcing term in the momentum equation. Another way to in- 
ject energy into the system is by adding an irregular external 
electromotive force Sc^t(x^€) to the system. The induction 
equation (^) is then replaced by 

9A 

— ^ u X B + r/V^A + 5oxt(a;, <). (54) 

This is in some ways reminiscent of turbulence driven by 
magnetic instabilities (e.g. the magnetorotational instability, 
buoyancy or kink instabilities). The direct injection of energy 
into the induction equation was also considered by Pouquet 
et al. (1976), for example, who assumed the forcing to be he- 
lical. The results of such a run are shown in Fig. Note 



4.1. Magnetic helicity flux from global mean-field models 

The only information available to date about the magnetic he- 
licity of the sun is from surface magnetic fields, and these data 
are incomplete. Vector magnetograms of active regions show 
negative (positive) current helicity on the northern (south- 
ern) hemisphere (Seehafer 1990, Pevtsov, Canfield, & Met- 
calf 1995, Bao et al. 1999, Pevtsov & Latushko 2000). From 
local measurements one can only obtain the current helicity, 
so nothing can be concluded about magnetic helicity. Typi- 
cally, however, current and magnetic helicities agree in sign, 
but the current helicity is stronger at small scales. 

Berger & Ruzmaikin (2000) have estimated the flux of 
magnetic helicity from the solar surface. They discussed ol- 
effect and differential rotation as the main agents facilitat- 
ing the loss of magnetic helicity. Their results indicate that 
the flux of magnetic helicity due to rotation and the actual 
magnetic field obtained from solar magnetograms is nega- 
tive (positive) in the northern (southern) hemisphere. Thus, 
the sign agrees with that of the current helicity obtained us- 
ing vector magnetograms. Finally, Chae (2000) estimated the 
magnetic helicity flux based on counting the crossings of 
pairs of flux tubes. Combined with the assumption that two 
nearly aligned flux tubes are nearly parallel, rather than an- 
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tiparallel, his results again suggest that the magnetic helicity 
is negative (positive) in the northern (southern) hemisphere. 

What is curious here is that even though different scales 
are involved (small and intermediate scales in active regions 
and coronal mass ejections, and large scales in the differential 
rotation), there is as yet no evidence that the helicity spectrum 
of the sun reverses sign at some scale. 

We now discuss the magnetic helicity properties of mean- 
field ail dynamos in a sphere. There are different reasons 
why looking at mean-field dynamos can be interesting in con- 
nection with helical turbulence. Firstly, magnetic helicity in- 
volves one wavenumber factor less than the magnetic energy 
and tends therefore to be dominated by the large scale field. 
Magnetic helicity should therefore be describable by mean- 
field theory. Secondly, mean-field dynamos easily allow the 
spherical geometry to be taken into account. This is mainly 
because such models can be only two-dimensional. The re- 
sulting field geometry does in some ways resemble the solar 
field. We recall, however, that at present there is no satisfac- 
tory model that explains both the equatorward migration of 
magnetic activity as well as the approximate antiphase be- 
tween radial and azimuthal field {BrBcj, < during most 
of the cycle; see Stix 1976), and is still compatible with the 
observed differential rotation (dVl/dr > at low latitudes). 
The sign of differential rotation does not however affect the 
sign of the magnetic helicity. In fact, the velocity does never 
generate magnetic helicity, it contributes only to the flux of 
magnetic helicity; see Eqs (^) and (^. 



Table 2. Properties of mean-field dynamos for different signs 
of a and radial shear, il' = dil/dr. All of these models are 
problematic when applied to the sun. Case (iii), where ao > 
and n' < 0, satisfies two important constraints (butterfly di- 
agram with equatorward migration, antiphase between radial 
and azimuthal field), but it is incompatible with the helio- 
seismology data which indicate that at least at low latitudes 
n' > 0. 





n' <o 


n' >0 


ao < 


Case (i) 
helioseismology not OK 
butterfly not OK 
phase rel. OK 
Ht^ and Qn negative 


Case (ii) 
helioseismology OK 
butterfly OK 

phase rel. not OK 
Hti and Qn negative 


ao > 


Case (iii) 
helioseismology not OK 
butterfly OK 
phase rel. OK 

Hn and Qn positive 


Case (iv) 
helioseismology OK 

butterfly not OK 
phase rel. not OK 
Hn and Qn positive 



For the axisymmetric mean field, B, of an afl dynamo 
we calculate the gauge-invariant relative magnetic heUcity of 



Berger & Field (1984) for the northern hemisphere (indicated 
by N), 

f (A + Ap) • (B - Bp) dV, (55) 

where Bp = V x Ap is a potential reference field inside the 
northern hemisphere, whose radial component agrees with B 
(see Bellan 1999 for a review). This implies that A^f, = Ap^, 
on the boundary of the northern hemisphere (both on the 
outer surface and the equator). 

In spherical coordinates an axisymmetric magnetic field 
can be expressed in the form B = S^t^ + V x A^cf). The 
potential reference field is purely poloidal. Therefore, the dot 
product of the toroidal components in the integral (^) is sim- 
ply 

= / (A^ + Ap^)B^ dV. (56) 

For the dot products of the poloidal components of the in- 
tegral we express the magnetic vector potential in the form 
A = A^4> + V X C^4>, where • V x V x (C^4>) = 'B4,. 
We can then write the integrand for the poloidal component 

as 

(V X • V X (A^ - Ap0)0 = 'B^(A^ - Ap0) 

-V- [Ax (3^-3p0)<^]. (57) 

(Note that Cp^ = 0, as the potential reference field is purely 
poloidal.) The integral over the divergence term vanishes be- 
cause A^ = Ap^ on the boundary, so we have 

ff^P°') = [ B^{A^-Ap^) dV, (58) 

and hence 

H^ = 2 [ A^B^ dV. (59) 

The time evolution of can be obtained from the mean- 
field dynamo equation (e.g. Krause & Radler 1980), which 
gives 

diJN/dt = 25'- Qn, (60) 
where S — Jj^ £ • B dV^ was already defined in Eq. 
£ — aB — titHqJ is the e.m.f. of the mean field, and Qn 
is the integrated magnetic helicity flux out of the northern 
hemisphere (i.e. through the outer surface and the equator), 
with 

Qn = -2^[(5 + Ux B) X -dS, (61) 

where U = r sin 9 il(f) is the mean toroidal flow. It will 
be of interest to split Qn into the contributions from the 
outer surface integral, Qns> ™d the equatorial plane, Qsq, 
so Qn = Qns + QEq- 

In Table 1^ we summarise the key properties of typical aH. 
dynamo models for different signs of a and il' = dfl jdr (cf. 
Stix 1976). Here we have also indicated the sign of the mag- 
netic helicity flux that leaves from the northern hemisphere, 
Qn. In all the models presented here the sign of Qn agrees 
with that of R^, which also agrees with the sign of a. This is 
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because, averaged over one cycle, {dH-^ / dt) eye vanishes, so 
the sign of (QN)cyc rnust agree with that of (SVyc which, in 
turn, is determined by the sign of a; see Eq. (pO|). 

What remains unclear, however, is whether the negative 
current helicity observed in the northern hemisphere is in- 
dicative of small or large scales, and hence whether the sign 
of a is positive (case iv) or negative (case ii), respectively. 
Only in case (iv) there would be the possibility of shedding 
magnetic helicity of the opposite sign; see Sect. 4.3 below. In 
case (iii), Hf^ is also positive, but here ft' < 0, which does 
not apply to the sun. In case (ii), on the other hand, is 
negative, so shedding of negative magnetic helicity would act 
as a loss. 

Our goal is to estimate the amount of magnetic helicity 
that is to be expected for more realistic models of the solar 
dynamo. We also need to know which fraction of the mag- 
netic field takes part in the 11 -year cycle. In the following 
we drop the subscript N on iJ and M. For an oscillatory dy- 
namo, all three variables, H, M, and Qjouic vary in an oscil- 
latory fashion with a cycle frequency uj of magnetic energy 
(corresponding to 11 years for the sun - not 22 years), so 
we estimate |di7/d<| ^ ll!\H\ and Qjouic ^ loM, which, by 
arguments analogous to those used in Sect. ^T^ , leads to the 
inequality 

■'-skin 1 

(62) 

where ^kin = (2??/w)^''^ is the skin depth, now associated 
with the oscillation frequency oj. Thus, the maximum mag- 
netic helicity that can be generated and dissipated during 
one cycle is characterised by the length scale |i/ |/(2/ioM), 
which has to be less than the skin depth £skin- 

For r] we have to use the Spitzer resistivity which is 
proportional to T~'^/^ (T is temperature), so 77 varies be- 
tween 10"* cm^/ s at the base of the convection zone to about 
10^ cm^/ s near the surface layers and decreases again in the 
solar atmosphere. Using oj — 27r/(llyr) — 2 y. 10^^ s^^ 
for the relevant frequency at which H and M vary we have 
^skin ~ 10 km at the bottom of the convection zone and 
^skin ~ 300 km at the top. This needs to be compared with 
the value |iJ|/(2/ioA/) obtained from dynamo models. 

In order to see how the condition ( |62[ ) is met by mean- 
field models, we present the results from a typical aVl model 
of the solar dynamo. Although mean-field theory has been 
around for several decades, the helicity aspect has only re- 
cently attracted attention. In the proceedings of a recent meet- 
ing devoted specifically to this topic (Brown, Canfield, & 
Pevtsov 1999) helicity was discussed extensively also in the 
context of mean-field theory. However, the possibility of he- 
licity reversals at some length scale was not addressed at 
the time. We recall that in the Babcock-Leighton approach 
it is mainly the latitudinal differential rotation that enters. 
We also note that, although the latitudinal migration could 
be explained by radial differential rotation, meridional circu- 
lation is in principle able to drive meridional migration even 
when the sense of radial differential rotation would otherwise 
be wrong for driving meridional migration (Durney 1995, 
Choudhuri, Schiissler, & Dikpati 1995, Kiiker, Riidiger, & 
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Fig. 13. Butterfly diagrams of the toroidal field from a mean 
field model, compared with butterfly diagrams of the mag- 
netic helicity density near the surface and magnetic helicity 
flux due to differential rotation at the surface. Note that the 
local magnetic helicity density is most of the time positive 
(light grey) in the northern hemisphere and negative (dark 
grey) in the southern hemisphere. The magnetic helicity flux 
shows a similar variation, but there is a certain degree of can- 
cellation taking place at any moment in time. Cq = 0.175, 
Cn = ly- 10^ (lat). 



Schultz 2001). We therefore proceed with the simplest pos- 
sible model in spherical geometry. An example of the result- 
ing butterfly diagrams of toroidal field, magnetic helicity and 
magnetic helicity flux is shown in Fig. |l3|. 

The main parameters that can be changed in the model 
are ao (we assume that a = ag cos 0, where is colatitude), 
77x (=const), and the magnitude of the angular velocity, 
in the expression 



1-0.038 



0.004 



(63) 



which is a fit that resembles closely the surface angular ve- 
locity of the sun (Riidiger 1989). What matters for the ax- 
isymmetric dynamo is only the pole-equator difference, Af2, 
which is about a third of JIq. Our model is made dimension- 
less by putting i? = r/x = 1- The governing non-dimensional 
parameters are Cq, — aoR/rjT and Co = AflR^/riT, where 
Afi = finiax — f^min- This Icads to the observables Spoio, by 
which we mean the radial field at the pole, which is about 2 
gauss for the sun, and -Bboit. by which we mean the azimuthal 
field in the main sunspot belts. The solar value of -Bbcit is es- 
timated to be around 4 ... 8 kG, based on equipartition argu- 
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Table 3. Results from a mean-field dynamo model in a spherical shell with outer radius R and inner radius 0.7 R, with either 
purely latitudinal shear '(l^t)' or purely radial shear '(rad)'. The values of a given here (in units of tit/R) lead to dynamos 
that are only sUghtly supercritical. The cycle period, Tcyc, is given in non-dimensional form and would coiTespond to 22 yr 
for ijT = 2.4 X 10^^ cm^/ s (for the models with latitudinal shear), or 3.7 x 10^^ cm^/ s (for the models with radial shear). 
The range of -Bpoic/^bcit that is consistent with observations, (2 ... 5) x 10~^, is highlighted in bold face. Both A/n and 
vary with the dynamo cycle period; H^'^^ and (i/N) refer to the maximum and time averaged values of H^, respectively. 
(QNs)Tcyc is the outwards magnetic helicity flux through the northern surface, while ((3Eq)7cyc is the magnetic helicity flux 
through the equator into the southern hemisphere. 
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ments and on the total emergent flux during one cycle.^Thus, 
we expect -Bpoic/i?bcit to be in the range (2 ... 5) x 10"^. 

From the model results shown in Table |] we see that, 
once Bpoio/Bboit is in the range consistent with obser- 
vations, (highlighted in bold face), {Hj^s) /{2fioM-[.^R) is 
around (1 . . .3) x lO"'' for models with latitudinal shear 
using Eq. (|63|). Given that R = 700 Mm this means that 

/ {2fiQM^) « 70 . . . 200 km, which would be compatible 
with £skin only near the top of the convection zone. 

In units of B^^^^^R'^, the magnetic helicity flux through 
the northern hemisphere is about (2 ... 5) x 10"''. Using typ- 
ical parameters for the sun, i?bcit — 4 kG and R = 700 Mm, 
we have B^^i^R^ = 4 x 10^° Mx^. Thus, our models pre- 
dict magnetic helicity fluxes of about (1 ... 2) x lO''^ Mx'^ 
for a full 22 year magnetic cycle, or half of that for the 1 1 
year activity cycle. This is comparable to the values reported 
by Berger & Ruzmaikin (2000), but less that those given by 
DeVore (2000). Significant magnetic helicity fluxes through 
the equator are only found if angular velocity contours cross 
the equatorial plane, i.e. in the models with radial differential 
rotation. 

If the observed negative magnetic helicity flux on the 
northern hemisphere is associated with the large scale field, it 
would suggest that ao < 0. This corresponds to case (ii), 
where magnetic helicity flux would contribute as an addi- 
tional loss of large scale magnetic energy. On the other hand, 
if the observed negative magnetic helicity flux is associated 
with the small scale field, it could remove excess magnetic 
helicity of the opposite sign, which may act as a source of 
large scale magnetic energy; see Sect. 1-.3. This corresponds 
to case (iv) where ao > 0. 



^ For the sun the total emergent flux per cycle is about 10^'' Mx. 
Assuming that the toroidal field is distributed over a meridional sur- 
face of ~ (20-50) X 500 Mm^ ^ (1-3) x 10^° cm^, this amounts 
to an estimated toroidal field of 3 — 10 kG. 



4.2. Balance equations 

The models with imposed positive helical forcing at small 
length scale show clearly that there is a production of equal 
amounts of positive and negative magnetic helicity at small 
and large scales, respectively (or vice versa if the helicity of 
the forcing is negative). As time goes on, more and more 
positive and negative magnetic helicity builds up on small 
and large scales. Of course, magnetic helicity always implies 
magnetic energy as well (this is related to the realisability 
condition; see Sect. 2.2 and Moffatt 1978), and eventually the 
value of the small scale magnetic energy will have reached 
the value of the kinetic energy at that scale, so the dynamo 
reaches saturation. This happens first at small scales, and if 
that was where the dynamo stops, the resulting field would be 
governed by small scales (Kulsrud & Anderson 1992). How- 
ever, this is not what happens in the majority of the numerical 
experiments, and probably also not in real astrophysical bod- 
ies like the sun. 

Before we go on, we need to address the question of op- 
posite helicities in the two hemispheres of astrophysical dy- 
namos. For reasons of symmetry, the total magnetic helic- 
ity is always close to zero. It is then convenient to think of 
the two hemispheres as two systems coupled with each other 
through an open boundary (the equatorial plane). The equa- 
tion of magnetic helicity conservation has an additional flux 
term, Q, on the right hand side, so it reads [cf. Eq. (|])] 

— = -2rjC - Q, (64) 
at 

where H is the magnetic helicity, C the current helicity, rj the 
microscopic magnetic diffusivity, and Q the magnetic helicity 
flux integrated over the surface of the given subvolume (i.e. a 
hemisphere in this case). Q We are interested in the evolution 
of the large scale field B, and we refer to the magnetic helicity 
of this large scale or mean field as Hm- 



Q is not to be confused with the rate of Joule heating, Q j( 
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In the absence of boundaries we have iJ,n = / A • B dV^, 
where B = V x A and overbars denote a suitable aver- 
age that commutes with the curl operator and satisfies the 
Reynolds rules (in particular, the average of fluctuations van- 
ishes and the average of an average gives the same; see 
Krause & Radler 1980). With open boundaries, however, 
/ A • B dV^ is no longer gauge invariant. For the present 
discussion we only need the fact that it is possible to de- 
fine gauge-invariant versions of magnetic helicity; see Berger 
& Field (1984) for the three-dimensional case and BDOl (or 



Sect. 2.4) for the one-dimensional case that is relevant to hor- 



izontally averaged fields. 

Splitting the field into contributions from mean and fluc- 
tuating components we can write 

^ + -TT^ = -'^vCrn - 2rjCf - - Q{. (65) 
at at 

This equation only becomes useful if we can make appro- 
priate simplifications. BDOl considered the case where the 
small scale field has already reached saturation and the large 
scale magnetic energy continues to build up further until it 
too reaches saturation. 

If the large scale field is nearly fully helical we have 
Mm w ■^k^\Hai\, where fcm is the typical wavenumber of 
the mean field. Furthermore, if the small scale field is nearly 
fully helical, too, we have Mf — ^kf\Hf\, where fcf is the 
typical wavenumber of the fluctuating field, or the forcing 
wavenumber, as in BOl. Analogous relations link C to H. 
The relative ordering of the ratios of large and small scale 
magnetic energies relative to magnetic and current helicities 
is given by 

(66) 



kf M„ 



\Hf\ k 
Thus, if k„ 



1 Mi 
< fcf, 



ICfl 



as is usually the case when there is some 
degree of scale separation, and Afm w il/f, we have |-ffm| ^ 
|iJf I and |C,„| < |Cf |, so we are left with 



dt 



(67) 



or, using the above estimates between M, H and C, and as- 
suming, for definiteness, that H^^, C,n > and Hf,C{ < 
(as expected in the northern hemisphere of the sun), we have 

« +2?7fc,„fcfMf - 2fc„Q,„ - 2fc,„Qf. (68) 

at 

(The factor 2 in front of the Qs comes from the 1/2 in the 
definition of M.) We can apply this equation to the problem 
with open boundaries. 

However before doing this, we recall that in the absence 
of any surface fluxes of helicity (Q — 0) the assumption 
Afni = M{ does not apply near saturation. Instead, steadi- 
ness of H immediately requires C = Cm + C[ ~ 0, so 
jCml = |Cf|, and hence M,n = {kf/k^)Mf > Mf, and 
l-f^ml ^ \Hf\ for the final state. Before this final state is 
reached we have from Eq. ( |65| ) 

dAfm 



dt 



-2r/fcmfcfMf - 2r/fc;^M„ 



(69) 



which has the solution (BOl) 

-2i)fc,^(t-iB 



k{ 



1 



(70) 



This equation is referred to as the magnetic helicity con- 
straint, which is probably obeyed by all helicity-driven large 
scale dynamos. Note, however, that no mean-field theory has 
been invoked here. Equation ( [70| ) implies that full satura- 
tion is reached not earlier than after a resistive time scale, 
(r/fcm)~^. We emphasise that this is the resistive time scale 
connected to the large length scale 2Ti/k-^. We also empha- 
sise that this result is unaffected by turbulent or other (e.g. 
ambipolar; see BSOO) magnetic diffusion. 

Coming back to the problem with open boundaries we 
can carry out a similar analysis by relating Qm to the mean 
field. Based on dimensional arguments we may assume 

IQml = Afm/(fcmrm), (71) 

where Tm is some damping time for the large scale field. We 
can express this in terms of an effective (turbulent-like) mag- 
netic diffusion coefficient, r/m, for the large scale field, e.g. in 
the form 



2ryn 



(Tm^m) ^ (definition of ?7,n), 



so Eq. ( p8[ ) becomes 



dt 



2ij^kiM„ 



^kraQl 



which has the solution 

2r/fcfMf - 2Qf 



2Tjin kn 



1 - 



(72) 



(73) 



(74) 



We expect i]^ to scale with the turbulent magnetic diffusiv- 
ity, rjt, which describes the spreading of the mean magnetic 
field. This equation yields the time scale {T]^k'^)~^, which is 
indeed shorter than the resistive time scale, (77/0^)^^- How- 
ever, if Qf = in Eq. (Q), the energy of the final field 
relative to the small scale field is decreased by the factor 
{rjkf) / [rji^k^^) ^ 1. This was the basic result of BDOl where 
I Qf I was found to be much smaller than | Qm | . In other words, 
all the magnetic large scale helicity loss occurred at large 
scales, which is of course the scale where we wanted to gen- 
erate the large scale field. 

The above results have shown that magnetic helicity 
losses due to large scale magnetic fields are not really a sat- 
isfactory solution to the problem of generating strong large 
scale magnetic fields on fast enough time scales. Rather, one 
requires that there should be a loss of magnetic helicity due 
to small scale magnetic fields, i.e. a non zero Qf. We may 
assume that Qf and Qm have opposite sign, because both 
current and magnetic helicities also have different signs for 
the contributions from small and large scales. Looking back 
at Eq. ( |74| ) is it clear that the Qf term in the numerator can 
in principle enhance and even completely supersede the re- 
sistively limited 'driver', 2rjkfMf, provided Qf < on the 
northern hemisphere. 

Physically speaking, the role of small scale magnetic he- 
licity flux is to let the large scale field build up independently 
of any constraints related to the small scales. To prevent then 
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excessive build-up of magnetic energy at small (and interme- 
diate) scales, which diffusion was able to do (albeit slowly), 
we need to have another more direct mechanism that is inde- 
pendent of microscopic magnetic diffusion. In order to in- 
vestigate the problem of excess magnetic helicity at small 
scales we have performed a series of numerical experiments 
by Fourier filtering the field and removing in regular time in- 
tervals the magnetic energy at the forcing and smaller scales. 

4.3. Removing excess magnetic helicity 

According to our working hypothesis, the reason why the 
large scale field can only develop slowly is that the build- 
up of the large scale field with one sign of magnetic helicity 
requires the simultaneous build-up (and eventual removal) of 
small scale field with opposite magnetic helicity. Instead of 
waiting for ohmic field destruction one might expect that also 
the explicit removal of magnetic energy at small scales would 
be able to accomplish this goal. 

A more formal argument is this. We use again Eq. (^5|), 
but without ignoring the dHf / dt term, and make the assump- 
tion that large and small scale fields are fully helical, and 
ignore surface terms, because boundary conditions are as- 
sumed periodic. Thus, we have an additional dMf/dt term 
in Eq. (p9|), which then reads 



dMm dM{ 



dt 



k{ dt 



(75) 



After small scale magnetic energy is removed, the right hand 
side of Eq. ( [75| ) is out of balance. If we ignored the dMf/dt 
term, the large scale magnetic energy would decay propor- 
tional to — 2ryfc^-jAfm, because the 2riknik{M{ term is still 
negligible after having set Aff = 0. However, the M[ term 
begins to recover immediately which gives rise to a signif- 
icant contribution from the dAff /dt term. Both terms enter 
with the same sign. If the magnetic Reynolds number is large, 
i.e. if the recovery rate of A/f is faster than than the micro- 
scopic decay time, we have a net gain. This is indeed borne 
out by the simulations, where magnetic energy is removed 
above the wavenumber fc = 4 every At — 0.12A"^ time 
units, labelled (a) in Fig. and At = 0.4A"\ labelled (b). 
Here, A = 0.07urmsfcf is the kinematic growth rate of the 
rms field strength. Shorter time intervals do not give any fur- 
ther enhancement; the resulting curve for At = 0.04A^^ lies 
almost on top of the curve (a) that starts at Xt = 80. 

We may thus conclude that the preferential removal 
of small scale magnetic fields does indeed enhance large 
scale dynamo action. Earlier experiments, which are not re- 
ported here, where k[ was larger, and therefore the magnetic 
Reynolds number Urms/?y^f smaller, did not show such an 
enhancement. Note however that we still do not really know 
whether the preferential loss of such small scale field does 
operate in real astrophysical bodies. Simulations have so far 
been unsuccessful in demonstrating this. Those presented in 
BDOl show that most of the magnetic helicity is lost through 
the large scale field. A possible reason may be that the large 
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Fig. 14. The effect of removing small scale magnetic en- 
ergy in regular time intervals At on the evolution of the large 
scale field (solid lines). The dashed line gives the evolution 

2 

of (B ) for Run 3 of BOl (where no such energy removal 
was included) in units of B^^ = /io/3o(u^). The two solid 

2 

lines show the evolution of (B ) after restarting the simula- 
tion from Run 3 of BOl at Xt = 20 and Xt = 80. Time is 
scaled with the kinematic growth rate A. The curves labelled 
(a) give the result for At — 0.12A^^ and those labelled (b) 
for A< = 0.4A^^. The inset shows, for a short time interval, 
the sudden drop and subsequent recovery of the total (small 
and large scale) magnetic energy in regular time intervals. 



scale field geometry was not optimal and that in other settings 
large scale losses are reduced. 

The intermediate scales are probably those where most of 
the magnetic helicity loss from the sun occurs. This argument 
is based in particular on the sign of the flux of magnetic he- 
licity on the solar surface: it is found to be negative (positive) 
on the northern (southern) hemisphere (see §4.1 above). 



According to conventional mean-field theory the sign of 
the a-effect should be positive in the northern hemisphere, so 
the sign of large scale magnetic helicity and the flux thereof 
should be positive. A reversal of the sign of magnetic helic- 
ity can then be expected on some scale smaller than the large 
scale of the dynamo wave (100-300 Mm). Thus, the helic- 
ity flux occurs on the scale of active and bipolar regions as 
well as the scale of coronal mass ejections, which is about 
20-50 Mm. These scales may therefore well correspond to 
the intermediate scales referred to above. On the other hand, 
there is some uncertainty as to whether the dynamo alpha in 
the northern hemisphere of the sun might in fact be negative. 
If the dynamo a is negative in the northern hemisphere, the 
magnetic helicity loss on the sun would be of the same sign 
as that of the large scale dynamo-generated field, so its scale 
would still belong to the large scales. However, the solution 
to the helicity problem would then probably not be related 
to magnetic helicity flux, because then Qf and Qm, as well 
as iJm, would have the same sign, so Qf would constitute 
an additional loss term, which does hence not accelerate the 
dynamo (BDOl); see also Sect. ^. 
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Fig. 15. Three-dimensional visualisation of the magnetic field for a simulation where the box is four times longer in the 
y-direction. Vectors in the back half of the box show the orientation of the locally averaged field. Ribbons in the front half 
indicate magnetic field lines; their twist reflects the local torsality (or helicity) of the magnetic field and has been enhanced 
artificially by factor of 15. The vertical magnetic field component is indicated by bright (positive) or dark (negative) colours 
in the top and bottom face of the box. 



4.4. Large aspect ratio simulations 

In order to assess the role of intermediate scales, we now ex- 
tend the simulations of BDOl with open boundaries to do- 
mains that are larger in one horizontal coordinate direction 
compared to before. We have already mentioned in the intro- 
duction that the simulations of BDOl produced most of the . 
magnetic helicity at large scale, i.e. the same scale as that of 
the large scale field that we want to build up fast. A potential 
problem with these simulations could be a lack of sufficient 
scale separation. In particular, it is conceivable that interme- 
diate scales are not sufficiently well represented in these sim- 
ulations. The intermediate length scales can be expected to be 
where most of the magnetic helicity of the opposite sign re- 
sides ('opposite' is here meant to be relative to the large scale 
field). 

In order to study the possible role of the allowance of in- 
termediate length scales in the simulations we extend the sim- 
ulations of BDOl to the case of larger aspect ratio. The hope 
is that the large scale field will now develop in the coordinate 
direction that is longest, but that most of the magnetic helicity 
flux occurs still at a similar scale as before and that the sign 
of magnetic helicity and its flux are now reversed relative to 
the helicity at the scale of the mean field. In the simulations 
of BDOl a box with unit aspect ratio was considered, and 
one of the coordinate directions (the z-direction) was non- 
periodic and magnetic helicity flux through these boundaries 
occurred. We have now extended this study to the case where 
the box is four times larger in one of the other directions (the 
y-direction). A three-dimensional visualisation of the result- 



Table 4. Summary of the main properties of the runs with 
vertical field boundary condition. The parameter qj gives the 
fractional magnetic energy in the mean field relative to the 
total magnetic field. 
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ing field is shown in Fig. |T5|. The y component of the mag- 
netic field on the z-boundaries shows a clear large scale pat- 
tern with a length scale larger than the extent of the domain 
in the x direction. The evolution of the magnetic and kinetic 
energies of one such run is shown in Fig. 16 We see that 



kinetic and magnetic energies are approximately in equipar- 
tition. The large scale magnetic energy shows, as expected, 
variation in the y direction, i.e. the wave vector of the large 
scale field points in the long direction of the box. 

Different runs are compared in Fig. |l^. All the models 
with open boundaries have the property of developing rapidly 
a mean field that varies in the y-direction, which is the long 
direction. When the aspect ratio A is small, however, the final 
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Fig. 16. Evolution of kinetic and magnetic energies. The 
magnetic energy contained in the large scale field (averaged 
in the x and z directions, and denoted by Ky) is shown as a 
dash-dotted line. The dotted and dashed curves denote the z- 
and x-dependent components, respectively. 
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Fig. 17. Comparison of the evolution of the large scale mag- 
netic field component that varies in the y-direction (the direc- 
tion in which the box is longest) for runs with different aspect 
ratios. A, and different values of ?7_3 = Note the de- 

velopment of a stronger large scale field for A = 1 than for 
A = 4. 



field is not that which varies in the y-direction, but instead in 
the z-direction. This is seen from Table ^ where we compare 
the relative strengths of the mean field. 



(76) 



where (B) ^ is the mean field averaged in the two directions 
perpendicular to j (cf. BOl). In all the cases with A ~ 1, 
is larger than and qy. This is not the case when A — A: 
here q^ and q.^. are very small. Nevertheless, as the magnetic 



Reynolds number is increased, qy decreases (from 0.33 to 
0.19 when 77 is lowered from 10^"^ to 5 x 10^^). Also, com- 
paring qy in the A — A case with q^ is the A ~ 1 case, the 
mean field is clearly less strong. 

Thus, we may conclude that larger aspect ratios favour 
large scale magnetic field configurations that vary in the hori- 
zontal (y) direction instead of the vertical (z) direction. How- 
ever, at late times the saturation field strength is still decreas- 
ing with decreasing resistivity. It seems therefore worthwhile 
pursuing this work by studying still larger aspect ratios and 
yet different field geometries. 

4.5. Fractional small scale magnetic helicity 

In practice the magnetic helicity will never be near 100%; 
typical values are around 3-5% for convective and accretion 
disc turbulence (Brandenburg et al. 1995, 1996). This has 
prompted Maron & Blackman (2002) to consider the effects 
of a forcing that has only fractional helicity (i.e. the realis- 
ability condition for the velocity is not saturated). They found 
that there is a threshold in the degree of helicity above which 
(large-scale) dynamo action is possible. However, these re- 
sults were obtained at a resolution of only 64^ meshpoints, 
and it is therefore possible, that this threshold effect is really 
the result of an increase of the critical magnetic Reynolds 
number above which the large scale a-type dynamo operates, 
and that the model with that resolution has dropped below 
this critical value for large scale dynamo action. Below we 
will estimate whether the dynamo number was in fact large 
enough for this to happen. 

The effect of fractional helicity of the forcing on the evo- 
lution of the large scale magnetic energy can be assessed in 
terms of the magnetic helicity constraint, using a generalised 
form of Eq. (|69|). In the case of fractional magnetic helicity 
the small and large scale magnetic helicities will be less than 
expected based on the magnetic helicity constraint, so 



\Cf\/kf « efMf 
and 

km I Hni I — I Cm | / fcn 



(77) 



(78) 



where ef < 1 and em < 1 denote the degree of the resulting 
helicity on small and large scales, respectively. Then, instead 
of Eq. (|70[), we have 

e[k[ 



-Mf 



1-e 



-2jjfef„(t-t,) 



(79) 



which is just ef / em times the expression on the right hand 
side of Eq. (^^. Two things can happen: if the small scale 
field is only weakly helical, i.e. Cf ^ 1, but still Cm ~ 1, then 
the large scale magnetic energy will be decreased. On the 
other hand, if the large scale magnetic field is only weakly 
helical, i.e. em ^ 1, but et « 1, then the large scale mag- 
netic energy will actually increase relative to the prediction 
of Eq. (70). In practice, both effects could happen at the same 
time, in which case the change of the large scale magnetic 
energy will be less strong. 
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Fig. 18. Evolution of kinetic and magnetic energies in two 
runs with fractional helicity. On the left hand panel, fcf = 5 
and T] — 2 X 10^^, whilst on the right hand panel, kf — 27 
and = 2 X lO"*. 



In Fig. 18 we show the results of two runs with frac- 



tional helicity forcing. In both cases the magnetic energy 
exceeds the kinetic energy by fef/Zcm multiplied by an 'ef- 
ficiency factor', ef/em- In the first case we have fct/Zcm = 5 
and Cf/em — 0.3, giving super-equipartition by about 1.5. 
This is roughly in agreement with Fig. |l8[ In the second 
case we have ki/k^ = 27 and et/em — 0.1, giving super- 
equipartition by a factor of about 2.7. Although the run in 
the second panel of Fig. [l^ has not been run for a full mag- 
netic diffusion time (which is here ten times longer than in 
the first case), it seems clear that strong super-equipartition is 
still possible. 

In the framework of dynamo theory (Moffatt 1978, 
Krause & Radler 1980), the large scale field is described by 
Eq. (36). In a periodic domain with minimum wavenumber 
ki the condition for large scale dynamo action is that the dy- 
namo number, Ca — a/ (77x^1), exceeds unity. Standard es- 



timates are a = — ir(ci; • u) and rjt 



(e.g., Moffatt 



1978), where r is the correlation time. For strongly helical 
turbulence with forcing at a particular wavenumber kf, we 



have I (a; • u) I 

^ _ \a\/ki 



etfct(u^ 
|a|/(rytfci 



so 



effcf/tfci, 



(80) 



?7 + 77t 1 + i?m 

where 1=1 + R:^^ is a correction factor involving the mag- 
netic Reynolds number, defined here as — rjt/rj. For the 
Runs shown in Fig. [l^, and with the numbers given in the 
previous paragraph, the values of Ca are 1.4 and 2.7, respec- 
tively. Thus, the dynamo number is supercritical. This is con- 
sistent with the clear growth of the large scale field seen in 
the simulations (Fig. |l8|). Maron & Blackman (2002) found 
that for fcf/fci = 5 the critical value of et is around 0.5. This 
corresponds to Ca = 2.5/t w 2, which should still be su- 
percritical. However, given that only modest resolution was 
used, the estimate for a may have been too optimistic. 

4.6. Fractional large scale magnetic helicity 

We have seen that the large scale field strength will decrease 
if the small scale field is not fully helical. However, accord- 
ing to Eq. ([79I), the large scale magnetic field strength should 
increase if the large scale field were less helical. In order to 
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Fig. 19. Evolution of kinetic and magnetic energies in a run 
with perfectly conducting boundaries in the z-direction. Note 
that the magnetic energy evolves in stages, similar to Fig. ^ 

demonstrate that this is true, we show in Fig. |9| the results 
of a run where the forcing is still 100% helical, but the large 
scale field is not fully helical. The latter has been achieved by 
adopting perfectly conducting boundary conditions, in which 
case it is no longer possible to have fully helical Beltrami 
waves. These calculations are otherwise similar to those pre- 
sented in BDOl, where the effects of nonperiodic boundary 
conditions was considered. 

In Fig. |l^ we see another remarkable feature that was al- 
ready seen in Fig. ||: the large scale magnetic energy grows 
in stages. Indeed, looking at the actual magnetic field pattern 
one sees that the number of nodes in plots of the mean field 
decreases with time, just like in Fig. ^. 

4.7. Position of the secondary peak 

In the original Run 6 of BO 1 there was a field of intermediate 
scale that grew fastest at fc = 7. This wavenumber was found 
to be in agreement with fc,nax = |a|/(277T), which is what 
one expects for an dynamo. Here yyx = rj + rjt is, the 
sum of microscopic and turbulent magnetic diffusivity and a 
represents the a effect that is caused by the helical small scale 
velocity. In the simulation, the only parameter that can be 
changed to make |a|/(2?7T) vary is the microscopic magnetic 
Reynolds number. In order to check our prediction regarding 
fcmax, we reduce the value of 77 from r] = 10~^ (as in BOl) 
t0 77 = 2xlO^''. We find that now the fastest growing mode 
is at fcinax ~ 15, which is, on a logarithmic scale, almost 
indistinguishable from the forcing wavenumber, kf = 30; see 
Fig. At larger scales, k < 7, the exponential growth is 
equally fast in all modes, independent of the (microscopic) 
value of J]. 

The reason for decreasing scale separation with increas- 
ing value of i?m is relatively easy to see. According to the 
dispersion relation for the dynamo in an infinite domain 
(e.g., Moffatt 1978) the wavenumber of the fastest growing 
mode is 

fcmax = \a\/2T]T = |a|/2t77t = effcf/2t. (81) 
So, unless R^ is small, in which case l = 1 + R^ can be 
large, fcmax will be around fcf/2 giving a scale separation of 
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Fig. 20. Scale separation in runs with different resistivity 
(first and second panel). Higher resistivity results in larger 
scale separation (third panel). 



1:2. However, for small values of gf, the scale separation can 
be much larger 

In the high i?m limit, a scale separation of just a factor of 
2 is hardly visible during the kinematic evolution of the field 
at that scale. Once the dynamo is in the nonlinear regime, i.e. 
once the small scale field has saturated, a will be quenched 
and fcmax must decrease further. 

In order to see that the sign of the magnetic helicity is in- 
deed different at small and large scales we split the field into 
right and left handed parts, B — B+ + (see Sect. 2.2 and 
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Christensson, Hindmarsh, & Brandenburg 2001), and show 
in Fig. |T]the magnetic energy spectra of B+ and B^, 
and , respectively. Inverse transfer is only seen in the 
AfjT spectrum; the spectrum is dominated by the forc- 
ing scale. 



5. Dynamos with shear 



10 



Fig. 21. Decomposition of magnetic field into positive and 
negative helical parts. Note that the magnetic energy at the 
forcing scale has increased between t — 250 and t = 350. 
The magnetic energy at large scales is dominated by helicity 
with negative sign (dashed curves). 



5.1. Evidence for magnetic Reynolds number 
dependence 

Finally, we turn attention to the case of dynamos with shear 
The presence of shear allows for an additional induction ef- 
fect that is responsible for producing strong toroidal field 
from a poloidal (cross-stream) field component. According 
to a result from mean-field theory such dynamos can pos- 
sess oscillatory solutions (e.g. Parker 1979). In a recent study 
by Brandenburg, Bigazzi, & Subramanian (2001, hereafter 
BBSOl), where a velocity shear of the form Uy{x) oc sinx 
was applied to helically driven turbulence in a periodic box. 
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Fig. 22. Butterfly diagram of By, continuation from the run 
of BBSOl, but with r] = 2 x IQ-'^. Run (iii) of Table |. 

oscillatory solutions were confirmed also for non-mean field 
dynamos. The outstanding question, however, is whether the 
resistive time scale, which is now known to affect the satura- 
tion phase of the dynamo (BOl, BBSOl), also affects the cycle 
period. This question could not be answered conclusively in 
BBSOl, because only one run with one value of the magnetic 
Reynolds number was considered. 

In the meantime we have accumulated data from another 
run that had twice the value of the magnetic Reynolds num- 
ber The flow is driven by a forcing, which consist of two 
components: a term varying sinusoidally in the x-direction 
with wavenumber k = 1 driving the shear and a term consist- 
ing of random Beltrami waves with wavenumber fcf = 5. As 
in BBSOl, the ratio of the magnitudes of the two components 
of the forcing function is about 1 : 100, resulting in a shear ve- 
locity that is in the final state about 50 times larger than the 
poloidal rms velocity of the turbulence. This is large enough 
so that we can expect to be in the so-called 'ail regime' 
where oscillatory solutions are the rule. 

Corresponding butterfly diagrams of the toroidally aver- 
aged toroidal field components are shown in Fig. The two 
diagrams are for two different x-positions, x = =F7r and at 
X = 0, where the shear S — dul^^'' /dx takes on its nega- 
tive and positive extremum, respectively. The simulation was 
restarted from the run of BBSOl at the time t = 6 000 and run 
until t = 13 000 with rj reduced by a factor of 2.5 with respect 
to BBSOl . It turns out that the resulting magnetic field suffers 
a major disturbance after the magnetic diffusivity is reduced 
by just a factor of 2. More importantly, a clear pattern of a 
migratory dynamo wave is now almost absent. Instead, the 
toroidal flux pattern appears to just wobble up and down in 
the z direction. 

One possible explanation for this result is that the dy- 
namo wave was just very close to the threshold between 
oscillatory and nonoscillatory behaviour This is not very 
likely however: the estimates of BBSOl indicated that the 
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Fig. 23. Butterfly diagram of By, continuation from the run 
of BBSOl, but with ry — 10~^ (five times larger than in 
Fig. H). Run (i) of Table |. 

dynamo numbers based on shear, Cn = S/{riTkf), is be- 
tween 40 and 80, whilst the total dynamo number (V = 
CaCn) is between 10 and 20 (see BBSOl), and hence Ca ~ 
a/(?7Tfci) ~ 0.25. Thus, shear dominates strongly over the 
a-effect (Co/Cq is between 150 and 300), which is typical 
for ail-type behaviour (i.e. oscillations) rather than a^-type 
behaviour which would start when Co / Ca is below about 10 
(e.g. Roberts & Stix 1972). 

Another possibility is that an oscillatory and migratory 
dynamo is still possible, but it takes longer for the system 
to settle onto this solution. Unfortunately the present simula- 
tions have become prohibitively expensive in terms of com- 
puter time that this run cannot be prolonged by much fur- 
ther at this time. Instead, we have calculated a model with 
larger magnetic diffusivity and found cycle periods that are 
now somewhat shorter than before (700. ..800 compared to 
1000. ..2000 in BBS02), but in diffusive units they are sim- 
ilar (0.7.. .0.8 compared to 0.5. ..1.0 in BBS02); see Fig. || 
This suggests that the period in this oscillatory dynamo is still 
controlled by the microscopic magnetic diffusivity. A more 
detailed comparison of various parameters with the other run 
presented above and with that of BBSOl is given in Table |[ 

5.2. Can magnetic helicity ride with the dynamo wave? 

We now want to address the question whether magnetic he- 
licity can vary in space and time such that it is actually un- 
changed in the frame of the dynamo wave as it propagates. 
The background of these suggestions is as follows. Although 
the sign of the magnetic helicity in each hemisphere remains 
the same (and is opposite in the other hemisphere), its value 
still changes somewhat as the cycle proceeds. It is conceiv- 
able that even a small fractional variation may prove incom- 
patible with the Parker-type migratory dynamo operating on a 
dynamical time scale. One possibility worth checking is that 
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Table 5. Summary of the main properties of the three- 
dimensional simulations with shear. Here, rj/{cski) is the 
magnetic diffusivity in units of the sounds speed and the 
wavenumber of the domain, and Wcyc — ZTr/Tcyc is the cycle 
frequency. In Run (iii) there is no clear cycle visible. 



Run 



(i) 



(ii) 



(iii) 



»?/(csA;i) 

i{Bl)/{Bl)yr^ 

Mo(J-B)/(b') 

LOcyc/ivkl) 



10"^ 
5 

30 
600 

4 

20 
0.018 
0.11 
8.. .9 



5 X 10"" 
10 
80 
1200 
6 

30 
0.014 
0.06 
6... 12 



2 X 10" 

25 

200 
2500 

20 

60 
0.008 
0.014 
> 10? 



» 



0.0 0,5 



1.0 1.5 3.0 2.5 



3.0 



« 








3.0 



3.0 



Fig. 24. Butterfly diagram of By for the run of BBSOl with 
1] = 5x 10^^, compared with the corresponding diagrams of 
Tiy and for the same time interval. Note the presence of a 
'torsional' oscillation in Uy, i.e. a superposition of the driven 
shear flow (second panel), but the absence of any similar pat- 
tern in Uz (last panel). Run (ii) of Table || 



the magnetic helicity pattern could propagate together with 
the dynamo wave. In that case the magnetic helicity would 
be conserved in a Lagrangian frame moving with the wave, 
but not in an Eulerian frame. If this transport is due to ad- 
vection, we would then expect systematic sign changes of Uz 
in a space-time diagram, which does not seem to be the case 
however; see Fig. p4[ 

What does seem interesting, however, is the presence of 
what is usually (in the context of solar physics) refeiTed to 
as a 'torsional' oscillation (Howard & LaBonte 1980). The 



toroidal flow pattern clearly traces the toroidal magnetic field 
pattern. More recently Howe et al. (2000) found another 
shorter period of 1 .5 years which has been seen in helioseis- 
mological data. If this can indeed also be interpreted as the 
result of a similar magnetic field pattern (which is as yet un- 
detected) then this would suggest the presence of multiple pe- 
riods in the solar dynamo wave (cf. Covas, Tavakol, & Moss 
2001). 

In the discussion above we had mainly bulk motions in 
mind that would result in a transport of magnetic helicity. 
This does not need to be the case because magnetic helicity 
can also be transported along twisting magnetic field lines. 
This was used in the approach by Vishniac & Cho (2001), 
who proposed a dynamo effect based on a non-vanishing di- 
vergence of the magnetic helicity flux. If magnetic helicity 
were to ride with the dynamo wave, then this could perhaps 
coiTespond to the anticipated magnetic helicity flux. Unfor- 
tunately, in the parameter regime currently accessible to sim- 
ulations the effect proposed by Vishniac & Cho has not yet 
been confirmed (Arlt & Brandenburg 2001). We may there- 
fore conclude that transport or advection of magnetic helicity 
by meridional flows remains a possibility, but has not as yet 
been verified numerically. 

5.3. Helicity loss in the presence of surface shear 

The calculations of BDOl had the shortcoming that large 
scale shear was absent. Of course, shear does not produce 
net magnetic helicity, but it can lead to a spatial separation 
which could be particularly useful if there are surfaces. Large 
scale shear in the equatorial plane would wind up a poloidal 
field and hence would lead to magnetic helicity of opposite 
sign in the northern and southern hemispheres. In helically 
forced turbulence simulations, the various effects have only 
been studied in isolation: shear and helicity in BBSOl, sur- 
face shear but no helicity in Arlt & Brandenburg (2001), and 
different heUcities in the northern and southern hemispheres 
in Brandenburg (2001b). 

Surface shear may allow for the direct production and 
transport of magnetic helicity of different signs in the two 
hemispheres on a dynamical time scale. Consider a differ- 
entially rotating star permeated by an initially uniform mag- 
netic field parallel to the rotation axis. The differential rota- 
tion will wind up the magnetic field in different directions in 
the two hemispheres, where the field lines describe oppositely 
oriented screws. Such additional production and transport of 
magnetic helicity on a dynamical time scale may be impor- 
tant for the magnetic helicity problem. On the other hand, for 
this idea to be relevant we have to have a pre-existing large 
scale poloidal magnetic field, so this issue does not seem to 
address the problem of resistively limited field saturation. 

In order to find out whether large scale surface shear is 
important we have considered a model similar to Run C2 of 
Arlt & Brandenburg (2001), i.e. with finite shear in a layer 
of thickness 2tt and a quiescent (non-turbulent) halo with no 
shear outside. As before, the shear is given by Uy{x), while 



22 



Astron. Nachr./AN XXX (200X) X 



the boundary to the halo is at z = zbvr. However, in con- 
trast to the nonhelical calculations of Arlt & Brandenburg, 
we have now included a helical forcing that is negative (pos- 
itive) above (below) the midplane. The results are in many 
respects similar to the runs of BDOl with a halo: around the 
time of saturation there is a strong negative burst of small 
scale magnetic helicity together with positive small scale cur- 
rent helicity. The resulting magnetic helicity fluxes both out 
of the turbulent domain as well as upwards on the two bound- 
aries are fluctuating about zero. Thus, the addition of shear to 
the halo runs of BDOl seems to have only little effect on the 
evolution of helicity fluxes. 

The evolution of the magnetic and kinetic energies is 
shown in Fig. ^ After about 0.08 resistive times the dynamo 
has entered the saturation phase which is then completed af- 
ter 0.15 resistive times altogether. This is now much faster 
than in the case of periodic boundaries. Fig. ^ shows images 
of the mean toroidal magnetic field. Especially near the lower 
boundary of the disc surface (dash-dotted line) one sees the 
ejection of magnetic structures. 

As can be seen in Fig. |7[ the magnetic field structure is 
markedly different at x = ivr compared to x « 0, where 
it is stronger and pushed to the disc boundaries. This field 
expulsion to the boundaries could well be the result of usual 
field migration that is expected in the presence of shear and 
helicity. In terms of mean-field theory with positive a in the 
upper disc plane and negative below, we expect field propaga- 
tion away from (towards) the midplane when shear is positive 
(negative). This is indeed consistent with our model where 
shear is positive (negative) at a; = (a; = ±7r). 

6. Conclusions 

In an attempt to clarify the problem of resistive versus dy- 
namical time scales in models of the solar dynamo cycle 
we have tied up a number of loose ends that had been left 
open after the first exploratory simulations of BOl, BDOl, 
and BBSOl. There are two distinct issues that may or may 
not be controlled by the resistive time: the time it takes for 
the large scale field to reach saturation and, once saturation 
is reached, the time scale on which the the large scale field 
can undergo cyclic variations as seen in the sun. 

Looking back at Fig. |l] in the introduction we can still ex- 
pect various approaches to be successful, although many oth- 
ers now seem to have been eliminated. One possibility was 
that the dynamo may still operate on a fast time scale if, like 
in the kinematic case, the wavelength of the dynamo wave is 
shorter than the extent of the system, and possibly only some- 
what larger than the scale of the forcing. The main problem 
would then be to keep the process fast and to prevent the dy- 
namo from developing scales as large as the box size. How 
this may happen in reality is not clear. It may well be con- 
nected to the geometry of astrophysical dynamos which have 
spherical or disc geometry and are not box-like. This corre- 
sponds to Alternative A in Fig. |l]. The argument against this 
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Fig. 25. Semi-logarithmic and double-linear plots showing 
the evolution of kinetic and magnetic energies (dotted and 
sohd lines), compared with the magnetic energy of the az- 
imuthally averaged mean field. All quantities are normalised 
by the kinetic energy of the cross-stream motions. E^^°^^ 
refers to the total kinetic energy including the shear The ra- 
tio of the kinematic growth rate to the resistive time at large 
scales is here \lr\k\ = 160. 
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Fig. 27. Butterfly diagram of the mean toroidal field of 
Fig. ^ at a; = ±7r and x = 0. The lines at z = ±7r mark 
the boundaries of the turbulent slab. 



possibility is that even at intermediate scales the evolution to- 
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Fig. 26. Loss of magnetic field through the upper and lower boundaries. Images of the mean toroidal field are shown at 
different times. Shear is positive (negative) at a; = (a; = ±tt). t] = 10~* and = 5 x 10^^. The lines at z = ivr mark the 
boundaries of the turbulent slab. 



wards larger scales (i.e. the motion of the secondary bump in 
the spectrum, cf. Figs |l^ and [ill) seems to occur on a resistive 
time scale. This can be seen from the reduction of the speed 
of the spectral bump, atrav, as the resistivity is reduced (Ta- 
ble |l]). So the forcing scale may have to be not much smaller 
than the scale of the dynamo wave. Other plausible mecha- 
nisms could be based on non-a effect dynamos. Examples 
are negative magnetic diffusivity effects (cf. Zheligovsky et 
al. 2001), the incoherent a-effect (Vishniac & Brandenburg 
1997), and the Vishniac-Cho effect (Vishniac & Cho 2001). 
These are all subsumed under Alternative B in Fig. |l| 

The remaining alternatives are all based on conventional 
a effect dynamos, so they have finite net magnetic helicity 
(although of opposite sign in the two hemispheres). In or- 
der for the dynamo to operate fast enough, resistive magnetic 
reconnection (with magnetic helicity destruction) has to be 
fast enough. This could only be the case near the surface 
where the microscopic magnetic diffusivity is sufficiently 
large. This corresponds to Alternative C in Fig. [l] 

The other possibility is that the dynamo may work 
through losses of magnetic helicity of opposite sign, possi- 
bly at intermediate scales, for example near the surface (e.g. 
in coronal mass ejections and active regions) or across the 
equator This is referred to as Alternative D. The simulations 



presented in Sect. 4.3 have shown that this is in principle pos- 



sible. The problem here is that, so far, none of these mecha- 
nisms have been seen to occur naturally in any of the sim- 
ulations investigated so far It is therefore important to move 
toward more realistic simulations, perhaps in spherical geom- 
etry, with strong vertical stratification which could contribute 
to keeping the large scale field within the convection zone via 
turbulent pumping, but allowing fields above a certain thresh- 
old to escape through the surface. 



The purpose of this review was to outline the present sta- 
tus of our understanding of nonlinear large scale dynamos. 
Things are developing rapidly with new simulations appear- 
ing every month. So far it has mainly been a process of elim- 
ination, because many different ideas have been around, and 
constantly more ideas are appearing. The least explored case 
is that of dynamos in realistic spherical geometries. 

Although a number of direct simulations of MHD tur- 
bulence in spherical geometry have been analysed (Oilman 
1983, Olatzmaier 1985, Olatzmaier & Roberts 1995, Drecker, 
Rudiger & Hollerbach 2000, Ishihara & Kida 2000) the issue 
of fast versus resistively limited growth has not yet been in- 
vestigated. This is indeed not an easy task. In order to distin- 
guish between the two possibilities, one would need to verify 
explicitly whether a large scale field is generated on resis- 
tive times, or on a time scale that is independent of and much 
shorter than the resistive time scale. In addition, if the models 
show cyclic variability (which is desirable for solar dynamo 
models), the cycle length needs to be (asymptotically) inde- 
pendent of the resistive time if the dynamo is fast. One reason 
why such numerical experiments are difficult is that only at 
relatively large magnetic Reynolds numbers (which require 
high numerical resolution) the resistive time scale becomes 
sufficiently long so that it can clearly be distinguished from 
the dynamical time scale. A relevant dynamical time scale is 
the ratio A/ {rikf) which has to exceed a value of around 20- 
30 before one can see that resistively limited growth or cycle 
periods have occuiTed (see BOl). At the same time the resolu- 
tion must not be too high because otherwise one may not be 
able to run the simulations for long enough before any large 
scale field has occuiTed. The same is true of hyperresistivity 
which tends to make the large scale resistive time scale ex- 
tremely long so that nothing can be said about resistively lim- 
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ited growth. This explains why the issue of resistively Hmited 
growth is not yet well understood in more realistic geome- 
tries. Nevertheless, it would be worthwhile reanalysing data 
from recent, high-resolution dynamo simulations in the light 
of resistive limitations on the duration of the saturation phase 
and, if applicable, the cycle period. 
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